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SUMMARY 

The fluid/thermal environment in the cove between the wing and eleven 
surfaces has been a concern throughout the design and initial operational 
phases of the Space Shuttle. Experimental and analytical investigations 
provided some understanding of the environment in the wing-elevon cove. This 
paper describes an incompressible finite element analysis of flow through 
straight parallel walls and curved parallel walls to determine the effects of 
cove geometry on the fluid/thermal environment. Results from this analysis 
agree qualitatively with experimental data. The centerline gas temperatures 
and cold wall heating rates are virtually identical for the two cases 
indicating that the slight curvature has little effect on the overall thermal 
environment. 


INTRODUCTION 

The gap between the Space Shuttle wing and eleven surface is closed by 
seals at the eleven hinge line to prevent through flow of hot boundary layer 
gases. The environment within the resulting sealed cove, as well as the 
environment in the cove should the seal leak, has been an area of concern 
throughout the design and initial operational phases of the Space Shuttle 
because the hot boundary layer gas flowing into the cove has the potential of 
damaging the thermally unprotected inner wing and eleven structures. 
Experimental invest! gat ions^l^ provided some understanding of the 
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fluid/thermal environment in a wing-elevon cove. A simplified partial 
cross-sectional view of the 1/2 in. wide by 12 in. long wing-elevon cove used 
in these studies is shown in fig. 1. The cove wall formed by the wing was 
instrumented with pressure orifices and thermocouples which were used to 
determine the cold wall heating rates. The gas temperature distribution in 
the cove was measured with beaded thermocouples. A finite difference 2D 
Navier Stokes analysis of the external hypersonic and internal subsonic 
flow provided a qualitative understanding of the flow phenomena associated 
with the cove, but further analytical studies were needed to provide a better 
understanding of the flow inside the cove. 

In an effort to better understand the effects of cove geometry on the 
flow field and cold wall convective heat transfer to the cove surfaces, an 
incompressible finite element analysis of steady flow in the cove was 
performed by modeling the cove as both straight and curved parallel walls. 

This paper presents brief descriptions of the finite element model and code 
methodology, analytical results, and a comparison of the predicted heating 
rates with experimental data. 

DESCRIPTION OF ANALYSIS 

Governing Equations. - Flow in the cove was mathematically modeled with 
the conservation of mass, momentum (Navier-Stokes), and energy equations in 
two-dimensions for laminar, steady-state flow of a viscous fluid. The 
previous analysis(^) and the experimental data(^) indicated the Mach 
number in the cove was less than 0.3; therefore the flow was considered 
incompressible and viscous dissipation was neglected. The fluid viscosity and 
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thermal conductivity were temperature dependent. The resulting governing 
equations are a set of four partial differential equations in the four 
dependent variables: velocity components U and V, pressure, P, and 

temperature, T. 

Finite Element Methodology. - The finite element program NACH0s(3)was 
used to obtain a numerical solution to the governing equations. NACHOS was 
written to solve relatively small problems of free and low-speed forced 
convection. The program uses 8-noded subparametric and isoparametric 
quadrilateral elements. The flow and energy equations, coupled by the 
temperature dependent properties, are solved iteratively using Picard or 
Newton-Raphson algorithms. They are decoupled by solving each set separately, 
alternating between the two and updating the solution with each iteration. 
Models were generated and displayed using the built in features. Boundary 
conditions were applied in global coordinates. Results were displayed using 
NACHOS graphics and an added post-processor. 

Finite Element Model . - The entire flow region from thermocouple TCI to 
pressure orifice POl (fig. 1) was modeled. However, because of element 
limitations in NACHOS, the changes in cross-sectional area and geometry beyond 
thermocouple TC4 were not considered. The cove was first modeled as straight 
parallel walls and then, to determine the effects of curvature, it was modeled 
as curved parallel walls. The straight cove, using a symmetry boundary at the 
centerline, was modeled with 340 subparametric elements (17 x 20) resulting in 
over 3400 unknowns. The curved cove, inside/outside radius of curvature of 6 
in./6.5 in., was modeled with 400 isoparametric elements (20 x 20) resulting 
in over 3600 unknowns. Elements in both models were concentrated at the inlet 
and walls where gradients are steep. 
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The boundary conditions applied to both models were determined from 
experimental data(^) for a cove leak area of 12.5 percent of the entrance 
area and the free stream conditions noted on fig. 1. The inlet conditions 
were a specified uniform velocity of 258 ft/s, calculated from the data, and a 
uniform temperature of 1940®R, measured by TCI. The exit conditions were a 
uniform pressure of 73 psf, measured by POl. The wall temperature (Ty^) was 
specified to be a constant 550®R. 



T, = 3350“R 

Fig. 1. - Partial cross-section of experimental cove model. 

RESULTS AND DISCUSSION 

Selected velocity and temperature profiles from the straight sided model 
are shown in fig. 2. The velocity profiles, fig. 2a, show that the flow is 
developing throughout the cove. The slight overshoot in the profile near the 
inlet is resolvable by mesh refinement, but has an insignificant effect on the 
downstream velocity and temperature. The temperature profiles, fig. 2b, show 
that the centerline gas temperature decreases with distance downstream as the 
gas loses energy to the cold wall. This decrease causes the gradients at the 
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wall, and thus the cold wall heating rate to also decrease with distance 
downstream. The same profiles for the curved sided model are shown in fig. 

3. The velocity profiles, fig. 3a, show the flow at the inner wall 
accelerates and flow at the outer wall decelerates. This is due to curvature 
induced variations in the pressure gradients between the inner and outer 
walls. Further downstream, only centrifugal forces act, causing the profiles 
to skew toward the outer wall. The temperature profiles, driven by convective 
forces, behave similarly. Except at the exit where the model is inexact, the 
temperature gradients at the inner and outer walls are nearly the same and the 
heating rates vary by less than 1 percent. 

The predicted cold wall heating rate (Qy^) distributions are compared with 
experimental data in fig. 4. The straight and curved wall models give nearly 
the same results. The finite element analysis captures the general trend of 
the experimental data in the region where the cove walls are parallel, but is 
low and in some cases by a factor of 2. The scatter in the data in this 
region could be caused by the measurements being in the range of the error 
band of the instrumentation. The difference in the data and the analysis 
arises from uncertainties in the boundary conditions applied to the finite 
element model, since they were obtained from experimental data. Since the 
finite element model did not include the change in curvature and area past 60 
percent of the cove length, the analysis does not predict the increase in the 
downstream heating rate found in the experiment. 

Some comments should be made on the lessons learned in applying the 
finite element method to a realistic aerospace problem and on desirable 
features of future finite element codes. For user convenience, model 
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generation and display, which includes boundary conditions shown on the model, 
and graphical representation of results are absolute necessities. Also, 
boundary conditions should be applied in a local element coordinate system so 
that complex geometries can easily be modeled. Approximately 10,000 unknowns 
are needed to accurately model the wing-elevon cove. Other than computer 
storage limits, the code should not have a self-imposed element restriction. 
The implicit algorithm, which requires reassembly of the matrices for each 
iteration, resulted in large execution times (5000 CP seconds). Execution 
time for this type of problem could be reduced by using a more efficient 
algorithm, perhaps an explicit transient formulation. 

CONCLUDING REMARKS 

An incompressible finite element analysis was performed to determine the 
effects of geometry on the fluid/thermal environment in a wing-elevon cove. 

To determine the effects of curvature, the flow in the cove was first modeled 
as straight parallel flow, and then as curved parallel flow. The results were 
compared to experimental data from a geometry which the finite element models 
only approximated. The curvature in the section of the cove where the walls 
are parallel has little effect on the overall thermal environment in the cove; 
the heating rates are virtually the same for both straight and curved flow. 

The analysis gives results which agree qualitatively with experiment in this 
region. More detailed modeling is required at the cove exit, where the change 
in curvature and area could have a significant effect on the fluid/thermal 
environment. The analysis showed that finite element programs with 
significantly larger solution capacity and faster execution times are needed 
to solve realistic aerospace flow problems. The analysis has provided a 
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fundamental understanding of the fluid/thermal environment in a portion of the 
wing-elevon cove, but further studies are needed for a detailed definition of 
the flow in the entire cove. 
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(b) Temperature Profiles 


Fig. 2. - Profiles for straight cove model. 



(a) Velocity Profiles 



Fig. 3. - Profiles for curved cove model. 



Fig. 4. - Cold wall heating rate distribution. 
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